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Abstract: If we find a light Higgs boson at the LHC, there should be many observable 
channels which we can exploit to measure the relevant parameters in the Higgs sector. 
We use the SFitter framework to map these measurements on the parameter space of a 
general weak-scale effective theory with a light Higgs state of mass 120 GeV. Our analysis 
benefits from the parameter determination tools and the error treatment used in new- 
physics searches, to study individual parameters and their error bars as well as parameter 
correlations. 
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1. Higgs physics at the LHC 

Over the coming years, the main goal of LHC physics is to understand the breaking of 
the electroweak gauge symmetry. This symmetry we usually assume to be spontaneously 
broken, i.e. hidden by the presence of a non-symmetric vacuum. Such a vacuum can be 
induced through a gauge invariant Higgs potential involving powers of H^'H, where H 
is a SU{2) Higgs doublet [||, |2|. One of the degrees of freedom of this doublet becomes a 
physical Higgs scalar, while the other three are absorbed as longitudinal degrees of freedom 
into the massive W, Z gauge bosons. In this model all particle masses are proportional to 
the vacuum expectation value v = 246 GeV, and the Higgs couplings to other particles as 
well as the Higgs itself are proportional to that particle's mass. 

Gauge boson masses mwz arise from dimension-four WWHH and ZZHH terms. As 
the pseudo-Goldstone modes of the Higgs doublet are parts of the VF, Z fields there is no 
additional free parameter in this construction. Both masses m\\r z are given by v times the 
gauge couplings g times a rotation by the electroweak mixing angle. This relation between 
the (measured) gauge boson masses and v allows us to determine v without observing the 
Higgs scalar. 

Fermion masses arise from the dimension-four operators H"^^. This operator can be 
accompanied by any real dimensionless number, the Yukawa coupling yj. Once the Higgs 
field develops a vacuum expectation value the combination UfXv becomes the fermion mass. 
Note that in contrast to the gauge sector we have not yet measured any Higgs couplings 
to fermions, with the possible exception of yt in the electroweak precision analysis |^. 
The link between the known fermion masses and the largely unknown Yukawa couplings 
can be offset for example by non-standard Higgs sectors. 

A crucial question at the LHC will be: if we observe a light scalar particle, how can 
we tell that this state is the physical Higgs boson and not something else or a mixture 
between a Higgs boson and another new state. A good example for a state which looks 
like a Higgs scalar and mixes with a Higgs scalar is the radion in theories with large 
warped extra dimensions ||5|. A more conservative modification of the Higgs sector in the 
Standard Model could be an extension, for example by another Higgs doublet. A theory 
which requires such an extended Higgs sector is supersymmetry, where we find two CP-even 
scalar Higgs bosons which mix, leading to a light and a heavy mass eigenstate 0]. In the 
decoupling limit the lighter of these mass eigenstates becomes the Standard Model Higgs 
boson, but for finite supersymmetric and heavy-Higgs mass scales its properties deviate 
from the Standard Model case. 

Our strategy for studying the Higgs sector is analogous to our studies of new-physics 
parameter spaces Q: there are many pieces of information on the Higgs sector which are 
not continuous. These include the spin of the Higgs boson, its CP quantum numbers, or 
the structure of a CP-even WWH vertex [^. In this paper we study the continuous Higgs 
parameters, i.e. we assume that these quantum numbers are known, or in other words, we 
assume that we know the operator basis we need to include in the effective Lagrangian 
describing Higgs physics at the LHC. However, for each of these operators we have to 
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determine its strength, to then compare the result with the Standard Model or alternative 
weak-scale theories. 

Such an analysis will test all accessible Higgs couplings to known particles for a given 
Higgs mass |ll|, which we assume to be 120 GeV. This will first of all include the 



gauge boson couplings, where as mentioned above we have little room to alter the relations 
of the gauge boson masses and the gauge couplings. In addition, we will measure the 
Yukawa couplings and compare them to the known fermion masses. The observable Higgs 
channels at the LHC involve gluon fusion as well as weak boson fusion combined with a 
large number of branching ratios. Such an analysis will be statistically challenging. Just 
as the high-dimensional parameter space analyses studying new physics in the LHC era, 
the analysis needs to focus on the proper treatment of all errors |8|, [T^. 

Some Higgs couplings will most likely only be measured at a later stage, involving 
integrated luminosities of at least 300 fb~^. One of them is the only second-generation 



Higgs decay we might observe at the LHC, the decay into muons |14, 15]. Owed to the 
outstanding m^^j resolution in particular of CMS the weak-boson-fusion channel might 
come close to a 5a signal using modern analysis tools. However, because of the lack of 
experimental studies and due to this small event rate we neglect all light-flavor Yukawa 
couplings in this analysis. 

Extracting the one-loop Higgs coupling to is similar to the two-photon channel 
in suffering from a reduced branching ratio, plus an additional leptonic branching ratio 
of the Z boson. We do not include it in our analysis, because of its small event rate for 
low-luminosity running. 

Involving physics beyond the Standard Model a light Higgs can decay to invisible 
particles. The issue in searching for this decay is triggering the events, which should be 
possible for weak boson fusion production ||l^. Then, we can extract an invisible Higgs 
signal from an invisible Z decay using kinematic distributions and a careful analysis of 
signal- free control regions. Note that because of the limited energy it is unlikely that such 
a signal would be mimicked by the production of for example two neutralinos or charginos, 
because phase space essentially only allows for the production of one particle in this channel 
at the LHC |17|. 



Last but not least, the arguably hardest-to-measure Higgs coupling can be regarded 
as the most interesting: if the electroweak symmetry is really spontaneously broken, this 
means there exists a Higgs potential. Such a potential can only exist if there is a Higgs 
self coupling. While the WlWlH coupling is induced by such a self coupling, it would 
be most desirable if we could measure the explicit self coupling Xhhh at the LHC |T8| . 
Finally, a measurement of the quartic self coupling Xhhhh could confirm the actual form 
of the Higgs potential |l^, but it is unlikely that we will measure this last and final Higgs 



parameter at any collider currently discussed, including the VLHC or CLIC [20]. We skip 
such high-luminosity channels which should not feed back into our leading parameter set. 
The analysis presented in this paper can easily be extended. 

The final assumption we make for the structure of the Higgs sector is CP symmetry. 
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While in many new-physics models, like super symmetry, it is actually hard to avoid ad- 
ditional CP violation we currently have no experimental hints for such complex phases in 
the weak-scale Lagrangian [^]. In other words, the apparent absence of additional CP 
violation from the TeV scale could be considered one of the worst problems of TeV-scale 
model building; it serves as a benefit for our Higgs parameter analysis. 

Any parameter study at the LHC or elsewhere needs to focus on the error analysis. 
First of all, this means that we have to include statistical and systematic errors including 
a full correlation matrix as well as theory errors. The latter are numerically challenging, 
because they are not Gaussian. From these errors we need to construct a likelihood map 
over the entire parameter space, for example using Markov Chains or simulated annealing. 
The techniques we employ to scan the high-dimensional parameter space can have a large 
impact on the numerical efficiency, but they should be equivalent in the final result. One 
major technical difference between these methods is the study of secondary maxima in the 
likelihood map. 

An exclusive likelihood map is well suited for example to study correlations between 
different model parameters, i.e. different Higgs couplings, as well as local properties of 
the likelihood map. However, there are good reasons to ask more specific questions. The 
most obvious are about the probability or likelihood distribution of one Higgs coupling 
{e.g. gwwn) or the correlation of two of them [e.g. gwwH vs grrH)- The answer to this 
question depends on how we ask this questions, and what methodology we want to employ. 

Bayesian probability distributions and frequentist profile likelihoods are two ways to 
study a parameter space, where some model parameters might be very well determined, 
others heavily correlated, and even others basically unconstrained. Both of them rely on 
an exclusive likelihood map as a starting point. Following the SFitter strategy we carefully 
compare the benefits and disadvantages of the frequentist and the Bayesian approaches for 
each model parameter. 

The approach of mapping measurements onto a high-dimensional parameter space for 
example realized in SFitter ||8| (or in the largely similar Fittino tool [^]) are completely 
general: model parameters as well as measurements are included in the form of model 
and measurement files and can simply be replaced. SFitter serves as a general tool to 
map typically up to 20-dimensional highly complex parameter spaces onto a large sample 
of highly correlated measurements of different quality including a proper error treatment. 
While this approach is obviously mandatory for high-precision analyses, it is also best 
suited to extract the maximum information from any kind of correlated measurements. 

2. LHC measurements 
2.1 Production modes 

At the LHC we rely on two main production processes for a Standard-Model type Higgs 
boson: through a top loop the Higgs couples to two gluons, which leads to cross sections 
of the order of 10 — 100 pb (37 pb for a 120 GeV Higgs). The production cross section is 
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known at next-to-leading order including all top-mass effects |22, 23|. At NNLO the gluon- 
fusion production rate is known for an effective ggH coupling [^], where the description 
of gluon radiation off the top loop might not reproduce all distributions correctly. Due to 
unknown higher orders we assume a theory uncertainty of 13% on the production rate after 
acceptance and background rejection cuts |^5|. This error includes the uncertainty on a 
jet veto for H — > WW decays |26]. For a realistic analysis including actual LHC data we 
would have to check in detail that cuts or analysis strategies do not lead to large additional 
errors, for example by inducing large logarithms |27, p^]. 



The second largest production cross section is mediated by weak boson fusion, involving 
two forward jets with transverse momentum around miy/2 and typical cross sections of 
1 — 10 pb (4.5 pb for a 120 GeV Higgs). From a QCD perspective this production channel is 
particularly clean, because there is no color flow between the two DIS sides of the process. 
This simple fact we can make use of by vetoing central jets between the two forward tagging 
jets The total rate as well as the key distributions for this signal are known to next- 



to-leading order [^], including (as it turns out: negligible) interference effects |31|. An 
interesting aspect is that due to the suppression of higher-order QCD effects, strong and 
electroweak NLO corrections are of the same size but opposite sign. For this production 
process we estimate the theory error to range around 7% |3C]. 

For weak boson fusion there remains a caveat: Higgs rate measurements include the 
survival probability of a central minijet veto [p9| ]. It is one of the reasons why weak boson 
fusion signatures have such a spectacular signal-over-background ratio. While we will be 
able to measure this veto probability in Z-|-jets production [^] for the signal we need to 
predict it [^] — including the underlying event. While surely a minijet veto will work in 
LHC analyses, in Section 5.4 we will check what happens to the Higgs sector analysis with 
an increased error. 

The production channels of gluon fusion and weak boson fusion are not independent. 
Higgs boson production with two jets will have contributions from both gluon fusion and 
weak boson fusion. With appropriate cuts, we can enhance one channel over the other 
but in our rate analysis we nevertheless need to include both. We do not introduce an 
additional error for this separation, assuming it to be already covered by the theory errors 
and the systematic errors. 



A third production mechanism is the radiation of a Higgs boson off top pairs. This 
signature, in spite of its sizable rate of up to 600 fb (450 fb for a 120 GeV Higgs), suffers 
from its high-multiplicity final state and hence large combinatorial errors. Essentially, 
any Higgs decay product will also be present in the top decays. Its rate is known to 
next-to- leading order ||3^, with a remaining theory error around 13%. 

Last, but not least, there might be hope to see a Higgs boson produced in association 
with aW or Z gauge boson. Its reasonable rate of 0.3 — 3 pb (2.2 pb for a 120 GeV Higgs) 
suffers from the additional leptonic branching ratio of the gauge bosons. Still, in particular 
in the context of modern jet-structure driven searches this production channel could well 
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turn out useful if we combine it with a Higgs decay with a large branching ratio. For this 
channel we assume a theory error of 7%. 



Throughout this analysis we will skip the analysis of Higgs pair production [36|. In 
the Standard Model this channel has a one-to-one correlation with the Higgs self cou- 
pling 1 18] and is strongly suppressed. We note, however, that in models with additional 
Higgs resonances the pair production of light Higgs bosons might well be visible |37]. 



2.2 Signatures 

The challenge for Higgs searches at the LHC lies in triggering and the extraction of the 
signal from backgrounds. In our analysis we use the following list of most promising 
signatures for models not too different from the Standard Model [38|. This list should not 



be considered the final work for the LHC, but simply reflects our understanding at the 
current stage of Monte-Carlo studies: 

H ^ bb — while for small Higgs masses this is the dominant decay channel, with 
branching ratios up to 90%, it is particularly hard to extract from QCD backgrounds. 
Recent CMS and Atlas studies show that combinatorial backgrounds make it very 
hard to observe this decay in tiH production, where it might never reach the 5a 
level pg|]. The best chance for this decay channel will probably be the WH/ZH 
production mode [^] with a subjet analysis |41], currently under study by both 
collaborations. 

H — > WW — even though two leptonic W decays only allow us to reconstruct a 
transverse mass of the WW pair this channel is arguably the most powerful search 
channel for a light Higgs boson. In the gluon fusion process, where the most important 
observables to cut on are correlations between the two leptons this decay signature 
can be extracted using the opening angle of the two leptons (being small if the two 
W bosons come from a resonant scalar) |2^. In weak boson fusion, this decay can 
be used even if one of the Ws is far off-shell [42| with a spectacular S/B ~ 0{1). 



Combined with the ttH production mode this decay might turn out useful for on-shell 
decays |^^, so we include all three production mechanisms. 



H — > ZZ — due to its spectacular fully reconstructable four-lepton final state this 
signature is usually referred to as the 'golden channel'. It works in combination with 



the gluon fusion [p4| as well as with the weak boson fusion [45| production process, 
but is somewhat statistically limited to larger Higgs masses. 

H ^ TT — for this channel to be useful we need to reconstruct the invariant mass of 
the two taus. It is well known how to do this in the collinear approximation, where 



the taus are fast-moving [46|. Unfortunately, because of the steep gluon densities, 
gluon fusion produces Higgs bosons too close to threshold, which means we cannot 
use this production mechanism. In contrast, in weak boson fusion the Higgs recoils 
against the two tagging jets, so the taus are boosted. That way we can even measure 
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the Higgs mass to ~ 5 GeV [47|. Further improvement we expect from an improved 



statistical distinction of the signal from the Z — > rr background |48|. For the tiH 
production channel with a subsequent decay to taus we are not aware of experimental 
studies, so we leave it out for now 



^ 77 — the decay to two photons has the advantage of being the only fully 
reconstructable channel for a very light Higgs boson, as preferred by electroweak 
precision data. The 77 mass resolution of 2 GeV or better compensates for the 
small branching ratio and makes this decay channel a promising candidate for an 
observation in gluon fusion production [^^, weak boson fusion production |51] and 



even ttH associated production [^2|. Moreover, this channels allows us to measure 
the mass of a light Higgs boson with a precision 0(100 MeV), depending on the Higgs 
mass. 

H invisible — the best chance to observe a Higgs decay to invisible particles |^] is 
in combination with weak boson fusion, namely looking for two tagging jets recoiling 
against missing energy [^]. Of course, the determination of the Higgs branching 
ratio to invisible new states relies on our understanding of the production rate and 
therefore on the proper combination with all other Higgs signatures. Due to the lack 
of a full-simulation analysis, we do not include this signal as a measurement. On the 
other hand, because this decay channel can be of crucial importance for new-physics 
searches, we include it as an additional free parameter in part of our analysis. 

Because of the recent developments in the H ^ bb channels we briefly mention other 
strategies to see this Higgs decay at the LHC: first, there might be a faint chance to observe 
the bb decay directly in weak boson fusion. Producing a single Higgs boson might lead to 
a 2.9 — 5.2(7 signal with 600 fb~^ of triggered data |54]. Secondly, we can look for a Higgs 



boson produced in association with a W boson in weak boson fusion. This signature has 
been estimated to lead to a 4.4(T signal with 200 fb~^ of data including a functioning central 
jet veto [55]. The channel has been reanalysed recently leading to a significance S/^/B of 



about 1.8 |5^ for 100 fb ^ . Both of these analyses we do not include due to their limited 
potential during low luminosity running and a lack of full simulation. Last but not least, 



we can look for a hard photon and a Higgs in weak boson fusion |57]. For this channel 
we might expect a 1.2o" significance for an integrated luminosity of 30 fb^^ with a signal- 
to-background ratio around 1/20. The inclusion of a mini-jet veto, studied qualitatively 



in Ref. ||53l, would lead to an improvement by a factor 2 [57|. When discussing the tiH 
production channel with a Higgs decay to bottoms we should therefore keep in mind that its 
impact could be improved significantly after including production or similar channels 
in weak boson fusion. 



Note that including signal and background contributions to counting experiments 
makes it obvious that the rate measurements discussed above only rarely fulfill the ideal 
condition S ^ B. Instead, for gluon- fusion signatures we find that the number of signal 
events tends to be significantly smaller than the number of background events, with the 
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positive aspect that for example in the — > 77 channel a side bin analysis allows us to 
very precisely subtract the background events. For weak boson fusion signatures the ratio 
S/ B typically looks better, but on the down side the event numbers are comparably small, 
which in some channels forces us into the Poisson region. 

In our numerical analysis we use predictions for the production rates as described in 
Ref. ||5^. These we table the points and interpolate between them. For the branching ratios 
we use Hdecay |l59|l , modified to allow for arbitrary Higgs coupling values. To estimate the 
effect of cuts, we compute an efficiency factor from our rates and the values given in 
Ref. [S0|. Compared to this study we reduce the signal rate for ttH,H bb by 50%, in 
agreement with recent Atlas ||, ||] and CMS |6|] studies. For the WH/ZH, H ^ bb 
channel including a subjet analysis we use the rates given in the original paper and 
check for an effect of a degradation through detector effects in Section 



2.3 Errors 

In order to obtain reliable error estimates for the Higgs sector parameters, a proper treat- 
ment of experimental and theory errors depending on their origin is mandatory. We follow 
the CKMfitter prescription or Rfit scheme . The complete set of errors includes statisti- 
cal experimental errors, systematic experimental errors, and theory errors. The statistical 
experimental errors are treated as uncorrelated in the measured observables. In contrast, 
the systematic experimental errors for example from the luminosity or from tagging effi- 
ciencies are fully correlated. 

As efficiency factors are usually determined from other observed channels, they really 
behave like statistical errors and should therefore be Gaussian. The numerical values for 
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Table 1: Signatures included in our analysis for a Higgs mass of 120 GeV. The Standard Model 
event numbers for 30 fb^^ include cuts The factor after the background rates describes how 
many events are used to extrapolate into the signal region. The last two columns give the one-sigma 
experimental and theory error bars on the signal. 
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systematic errors we take from Ref. [S0| 



For convenience, we have included them in Table |2|. 



In contrast, theory errors reflecting unknown higher orders in perturbation theory 
should not be Gaussian but flat box-shaped within a certain range of validity of perturba- 
tion theory. In other words, the probability assigned to any measurement does not depend 
on its actual value, as long as it is within the interval covered by the theory error. There 
could be a tail attached to these theory-error distributions, but higher-order corrections 
are definitely not allowed to become arbitrarily large. Confronted with a perturbatively 
unstable observable one would instead have to rethink the perturbative description of the 
underlying theory. The numerical size of the theory errors we give in Table ^. The branch- 
ing ratios to cc and gg do not enter in any measurements but indirectly via the total width 
of the Higgs boson. To be consistent with our approach of flat theory errors the propagation 
of errors is linear and not quadratic. 

The error due to a mini-jet veto survival probability in background processes behaves 
like any other efficiency, i.e. Atlas and CMS measure it in signal-free phase space regions 
and assign a Gaussian error. However, for the signal rate we have to rely on theory to 
predict the survival probability. Hence, it is a fully correlated flat error which we can add 
to the error on the production cross section in weak boson fusion. 



Defining a scheme for flat theory errors includes their combination with the (Gaus- 
sian) experimental errors. A simple (Bayesian) convolution leads to the difference of two 
one-sided error functions with a clear maximum, so the convolution knows about the cen- 
tral value of theory prediction. A better solution is a distribution which is flat as long 
as the measured value is within the theoretically acceptable interval and then drops off 
like a Gaussian with the width of the experimental error, the Rfit profile likelihood con- 
struction p^]. The log-likelihood = — 21og£ given a set of measurements d and in the 
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Table 2: Systematic errors used in our analysis. Left: systematic errors applying to both signal 
and background. Reconstruction and tagging efficiencies are defined per particle, e.g. _ff 77 
has a 4% error on the photon reconstruction. Right: systematic background errors, either fully 
correlated or independent between channels. Tables same as Ref. Q. 
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presence of a general correlation matrix then C reads 
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'o 



Xd,i 



\di - di\ < a] 



(theo) 



dj- di + g, 

^(exp) 



(theo) 



di — di < —a. 



(theo) 



(2.1) 



di - d j - g- 

(exp) 



(theo) 



di-di > cr ■ 



(theo) 



cr- 



where di is the i-th data point predicted by the model parameters and di the actual measure- 
ment. Flat errors lead to a technical complication with hill-climbing algorithms. Functions 
describing box-shaped error distribution will have a discontinuities of higher derivatives. 
The Rfit scheme has a step in the second derivative which we need to accommodate. As 
a second complication the log-likelihood becomes constant in the central region, so some 
parameters vanish from the counting of degrees of freedom Q. If all experimental errors 
are Gaussian, (t(*^^p) of eq. {2A) is the convolution of the errors, i.e. the square-root of the 
sum of the squares of the experimental errors. The off-diagonal elements of the correlation 
matrix C are constructed by a generalization of eq. (2) of |p . 

The uncertainty on the mass of top and bottom quarks is another source of errors which 
we need to take into account. For these measurements we use a single experimental error 
each which comprises all different sources and is uncorrelated to any other measurement. 
Its numerical size is 1.0 and 0.07GeV for top and bottom quarks, respectively. 



For the experimental errors the systematic parts are inferred from large data samples, 
and it is safe to assume them as Gaussian. On the other hand, in some of the Higgs 
channels we will only observe a few events, so we need to use a Poisson distribution to 
model the behavior correctly. Therefore, we use Poisson statistics for all statistical errors. 
Our method of combining Gaussian and Poisson errors affecting the same measurement we 
describe in detail in Appendix Unless explicitly stated otherwise, SFitter uses smeared 
toy measurements to evaluate the errors. For each toy measurement we determine the 
best-fit value, which means the width of the distribution of the best-fit values gives the 
error on the model parameters. 
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Table 3: Theory errors used in our analysis for a 120 GeV Higgs. 
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3. General Higgs sector 



As sketched earlier, we allow for a general Higgs sector including couplings to all Standard- 
Model particles and an invisible final state. We assume that we observe a narrow 120 GeV 
Higgs candidate as a peak either in the 77 or ZZ invariant masses, in the collinearly 
reconstructed invariant rr mass, or in the transverse mass of WW pairs. This most 
notably means that we are searching for two-particle decays of the Higgs boson, which has 
an important consequence: if the Higgs boson were to couple strongly to light-flavor jets 
(which are swamped by QCD backgrounds) this would reduce the fraction of Higgs bosons 
decaying into visible channels, but it would also enhance the production cross section at 



the LHC. We discuss such a scenario as well as invisible Higgs decays in Section |4^. As 
mentioned in the introduction, we will not take into account Higgs couplings which might 
be visible in high-luminosity running in a single channel and/or are not expected to feed 
back into the determination of the leading Higgs parameters. This includes the Higgs- muon 
Yukawa coupling, the Higgs self coupling, or invisible Higgs decays. Such channels can be 
very useful to further understand the Higgs sector, but their results can be considered 
purely additional to our analysis. 

For the tree-level Standard Model Higgs couplings to any particle j we allow for a 
deviation 

am^ijfu (i + A,,/f) (3.1) 

where the Ajj// are independent of each other. If necessary redefining our Higgs field we 
can without loss of generality assume that one of the tree-level Higgs couplings is positive, 
i.e. gwWH > or A^yiy// > —1. Note that these couplings are of mass dimension three 
for gauge bosons and four for fermions, so they are part of the renormalizable Lagrangian. 
In the minimal one-doublet Higgs sector such corrections spoil gauge invariance. This can 
be cured for example by additional heavy Higgs states which we integrate out to define an 
effective one-Higgs-doublet Standard Model with free couplings. 

In addition, there are three relevant loop induced couplings in the Higgs sector: gggu-, Q'yyH 
and Q'yZH- Such couplings are sensitive to new particles in the spectrum, like for exam- 
ple supersymmetry or a chiral fourth generation. In the latter case the Tevatron limits 
on gg ^ H ^ WW strongly constrain the allowed Higgs masses in the presence of such 
a heavy generation which leads to an enhancement of the production rate roughly by a 
factor 9 |64]. Including a general new-physics scale A ||6^, ll^] such couplings arise in the 



effective Lagrangian as 



{GTiGnf.,^ + ^ F^'Ff.u (3.2) 



Note that in general new-physics models these couplings do not have to be perturbatively 
suppressed, as they are in the Standard Model or in the MSSM. In our weak-scale analysis 
we consider two different sources of higher-dimensional couplings: first, there are known 
Standard Model particles propagating in loops and inducing them, so any change in the 
tree- level couplings eq.(O) propagates into the one- loop effective couplings. Secondly, we 
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allow for manifestly dimension-five operators from new physics 

am a^jh (i + + A,,^) (3.3) 

where S^jh is the loop-induced coupling in the Standard Model, A^jf/ is the contribution 
from modified tree-level couplings to Standard-Model particles, and ^jjH is an additional 
dimension-five contribution, for example from new heavy states. 

The ansatz described above describes our general effective theory of the Higgs sector, 
allowing for deviation from the Standard Model values. Such deviations will be loop induced 
even in the Standard Model. If we allow for new physics in electroweak symmetry breaking 
such loop effects can be enhanced by large loop corrections, a strongly interacting sector 
mimicking a scalar Higgs boson, or an extended Higgs sector. We can illustrate this using 
9WWH- if we were to measure a smaller value than the Standard Model prediction this 
could point to additional Higgs fields sharing the burden of creating the correct gauge boson 
masses v'^ = Vi+V2 + - ■ ■ This kind of behavior we know from the type-H two-Higgs-doublet 
model of the MSSM where the couplings of the two Higgs scalars to the weak gauge bosons 
are split according to sin^(/? — a) + cos^(/3 — a) = 1. 

A too large value of gwwH, in contrast, could point towards the existence of an 
additional higher-dimensional operator. The way to check such a hypothesis would involve 
an analysis of angular correlations which can distinguish different tensor structures of the 
WWH coupling in weak boson fusion, independent of the decay channel Q. 

One basic assumption we have to make for the Higgs sector analysis is the treatment 
of the total width, which can be varied just like the couplings 

rtot = (1 + Ar) . (3.4) 

Allowing for a simultaneous scaling of all couplings and the total width we see that as long 
as we keep g'^/Ttot = C finite and constant, none of our {a ■ BR) measurements will be 
affected and the shift in the couplings will not be observable. Keeping this ratio and hence 
all measurements unaffected implies that for example a theory-motivated upper limit on 
any of the relevant couplings translates into an upper limit on the total Higgs width, and 
vice versa. 

Limiting ourselves to Higgs singlets and doublets we can derive such an upper limit 
on QwwH from the perturbative unitarization of WW WW scattering. As seen in the 
last paragraph, the unitarization of this process can be shared by several Higgs states, but 
since their couplings to gauge bosons enter twice per amplitude it leads to strictly a sum 
of the individual contributions. In other words, no individual WWHj coupling of any of 
the Higgs states j contributing can be larger than the value which alone unitarizes this 
process, i.e. g^^^j- This upper limit on gwwH can be translated into an upper limit on 
the total width. As mentioned in the last paragraph, this argument is limited to dimension- 
three WWH operators and for example breaks down in the limit of strong dimension-five 
interactions as shown in eq.( |3.3| ). 
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In a similar manner, we can obtain a lower limit on the Higgs width: obviously, the 
total width is at least as large as the sum of all observable partial widths. However, 
these partial widths are computed from the extracted couplings, and the couplings can 
be dialed to zero as long as the total width vanishes simultaneously. We can insert the 
scaling behavior Tvis oc 5^ into the expression for C, add an additional positive contribution 
Txig) > with an arbitrary g dependence and compute the small-coupling limit 

4 4 

C = lim ^= lim ^r-- — — — = . (3.5) 

However, by definition this ratio C has to be finite to protect our rate measurements, 
which means that we cannot actually dial down all partial widths simultaneously. Instead, 
we have to consider a quadratic equation in which indeed returns a lower limit on the 
couplings and hence the total width. It is obvious for any such argument that the knowledge 
of as many as possible partial widths and therefore the observation oi H ^ bb decays is 
crucial for studies like this. 



For a light Standard- Model Higgs boson, the difference between the observable partial 
width at the LHC and the Standard Model width is of the order of < 10%. Identifying 
the total width with the sum of the observed partial widths is practically equivalent to 
including an upper limit. Therefore, throughout this analysis we always simply identify 
the total width with the sum of the observed partial widths. Unless explicitly mentioned 
this includes all Standard Model particles listed in Section 0. An example scenario how 



this might mislead us we discuss in Section 4.3 



4. Higgs likelihood map 

Due to the number of possible observable channels for a light Higgs boson, the LHC should 
be able to not just discover a narrow scalar state which might or might not be the Standard- 
Model Higgs boson, but tell if such a state has anything to do with the Standard Model 
Higgs boson. This is a statistical statement, in which all errors and the resulting confidence 
level for a given hypothesis play a crucial role. Therefore, it is vital to control the mapping 
of the high-dimensional space of correlated measurements onto the high-dimensional Higgs 
parameter space. Just like for any other parameter study, the starting point of such a study 
is a completely exclusive likelihood map. 



Such a likelihood map is a large list or array of points in model space with their 
corresponding probability to agree with the given data. We can construct it using Markov 
chains where the points in the Markov chain represent the entire parameter space in their 
probability to agree with a data set. In principle, we could evaluate any kind of function 
on this Markov chain, similar to general mass scales or dark matter properties in analyses 
of physics beyond the Standard Model [37, 69]. In this analysis we are only interested 
in the probability to agree with data itself which in the case of supersymmetric parameters 
led us to consider weighted Markov chains 
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The problem with such an exclusive likelihood map is that it is (a) hard to graphically 
represent and (b) hard to interpret once we for example want to know if one individual 
coupling agrees with the Standard Model or does not. However, to reduce the dimension- 
ality of the parameter space we have to decide on a scheme, namely Bayesian probabilities 
or profile likelihoods. 

If we want to maintain the mathematical properties of a probability distribution and 
compare the height and size of different peaks we need to integrate over the unwanted 
directions, which means we need to define an integration measure in model space. Such a 
measure is always a free choice, a prior, and the resulting Bayesian probability distribution 
will depend on it. The model space integration also leads to volume effects, just like the 
usual example — is a water molecule more likely to be found in a highly concentrated 
spoon of water or in a huge cloud? Moreover, the integration will produce noise effects 
which can obscure narrow structures in the parameter space. 

Alternatively, we can construct a likelihood without introducing any measure, for ex- 
ample by simply picking the best point in the unwanted direction and identify it with the 
value in the reduced likelihood map. The result is referred to as a profile likelihood, has 
no noise effects and is designed to best represent the structures of the original likelihood. 
However, it does not allow us to compare the integrated size of different peaks of alternative 
solutions in the parameter space. 

In our new-physics parameter analyses [^] we have learned that we can extract the 
most information by showing both, a Bayesian probability as well as a profile likelihood, 
side by side, accepting that they do not have to give the same answers because they are 
not the same question. In this analysis we will find that correlations around a well-defined 
central fit value are more important than distinct equally good solutions, which implies 
that we will rely largely on profile likelihoods. 

Due to the probability distribution function which chooses their next point Markov 
chains are intrinsically Bayesian. Correspondingly, we find that profile likelihoods are 
harder to extract. To improve this behavior we can cool down the Markov chains to 
improve the resolution of local structures, provided we ensure that we do not miss structures 
in parameter space. If we slowly dial down the random number which decides whether a 
worse point is still accepted, we start scanning the parameter space similar to a standard 
Markov chain and then slowly concentrate on one structure. This cooling significantly 
improves the resolution of local structures around a peak and thereby yields a much better 
resolution for profile likelihoods. More details of this approach we give in Appendix 

4.1 Parameters and correlations 

Given the set of measurements described in Section ^ it is obvious that most of the Standard 
Model couplings should be accessible to a full analysis. Nevertheless, we start with a 
minimal set of Higgs sector parameters in which we only allow for tree-level couplings to 
all Standard Model particles. This implies that there are no new particles contributing to 
the effective ggH and 77// couplings. Since we compute the Higgs width as the sum of all 



-14- 




1 2 

^WWH 




1 2 3-5 -3 -1 1 3-5 -3 -1 1 

^WWH ^ttH ^bbH 



1 Q 1 2 3-5 -3 -1 13-5 -3 -1 1 

^WWH ^ttH ^bbH 







1 2 3-5 -3 -1 1 3-5 -3 -1 1 3 

^WWH ^ttH ^bbH 



10 12 3-5 -3 -1 13-5 -3 -1 1 3 

^WWH ^ttH ^bbH 



Figure 1: Profile likelihoods (left) and Bayesian probabilities (right) for the WWH, ttH, and bbH 
couplings. Not allowing for additional ggH or ^^H couplings we show results for 30 fb~^ and for 
300 fb~^ in the upper and lower rows. The Higgs mass is chosen as 120 GeV. All experimental 
and theory errors are included. Here and in all other Figures we assume the WWH coupling to be 
positive, i.e. Awwh > di- 



visible partial widths, a measurement of the bottom Yukawa constitutes the main fraction 
of the Higgs width. 

Based on the studies of weak boson fusion we limit our study to low-luminosity running 
and a conservative integrated luminosity of 30 fb~^. We can easily test how constraining 
this assumption is for our analysis: without including any effective higher-dimensional 
couplings we estimate the effect of an increase in luminosity from the standard 30 fb~^ to 
300 fb^^. Fig. m shows that the existence of correlations is independent of the integrated 
luminosity, i.e. these correlations are physical, and even the relative weights of alternative 
solutions hardly depend on the luminosity. Note that just like the entire analysis the error 
on the Higgs rates does not decrease like 1/y/S but that we always include backgrounds, 
which makes S + B the relevant number of events. 

Moving on, one striking feature in Fig. || is the alternative minimum for a negative top 
Yukawa coupling. In the comparison between profile likelihood and Bayesian probability we 
see that while the best-fit point clearly prefers the correct sign of the top Yukawa, volume 
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Figure 2: Correlations in profile likelihoods (left) and Bayesian probabilities (right), not allow- 
ing additional effective couplings. All experimental and theory errors included for low- luminosity 
running. 

effects wash out this distinction in the Bayesian probability. The only term which breaks 
this sign degeneracy is the loop-induced g^-yH, which only allows for a synchronized sign 
flip of guH and gwwH- The positive correlation between these two couplings shown in the 
two-dimensional panel of Fig. |l| arises on the one hand because of the strong destructive 
interference of the two loops and on the other hand due to a possible approximate scaling 
of all tree-level couplings compensating a shift in the less well measured bottom Yukawa 
coupling. A gwwH — gtbH correlation shows the same kind of positive correlations for the 
same two reasons. As the gbbH Yukawa coupling is small, a small shift in gwwH or guH 
compensates a major change in gtbH- However, such a large shift in the bottom Yukawa 
coupling is not completely out of the world — in the supersymmetric two-Higgs-doublet 
model the additional parameter tan /3 plays exactly this role. 

From the one-dimensional histograms in Fig. |l] we see that the LHC provides a very 
reasonable measurement of the top Yukawa coupling via the effective ggH coupling. The 
bottom Yukawa coupling is better constrained than in previous analyses using the tiH 



production channel, if we take the sub jet analysis at face value |41|. The measurement 
of gwwH is a little less precise than naively expected, given the weak-boson-fusion sig- 
natures including ample information on this coupling. Comparing the results for the two 
luminosities we see that the limiting factor for the WWH coupling is the statistics of the 
weak-boson-fusion channels. For the gluon-fusion production process combined with a de- 
cay H — > WW we know that the kind of off-shell decay appearing for light Higgs masses 
is essentially impossible to extract from the sizable backgrounds. 

Comparing the left and right set of panels, we see that in the Bayesian probabilities 
volume effects completely hide any kind of sign preference which might come out of the 
interference in the effective g'yyH and gggH- 

To highlight the importance of two-dimensional correlations, we show an extended set 
of them in Fig. |2|. In contrast to the correlations shown in Fig. |l| these are not induced by 
effective couplings, i.e. through quantum effects, but purely based on the (cr- BR) structure 
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our LHC measurements. 



First of all, there will be a general correlation between all couplings. Similar to scaling 
the total width as described in Section |^ we can shift the bottom Yukawa couplings which 
dominates the total Higgs width in the denominator of all rate measurements. We can 
compensate for such a shift by simultaneously scaling all relevant couplings, which gives a 
general tendency of a positive correlation of all couplings in this analysis. 

The correlation between gwwH and OttH is induced by the two weak-boson-fusion 
signatures, where from the same production process we can observe both decays at the 
LHC. This means, that the ratio gwwH / grrH is well constrained, without the possibility 
of compensating a shift of only one of the two decay couplings with another shift in the 
production cross section. This adds to the positive diagonal correlation gwwH oc grrH 
in the profile likelihood panel in Fig. ^. There will be a slight distortion due to the fact 
that we can keep the weak-boson-fusion — > rr rate constant with a negative correlation 

gWWH OC I /grrH- 

The correlation between gwwH and gzzH is less pronounced, mostly due to low statis- 
tics. In the weak-boson-fusion production process, the ZZ fusion channel typically accounts 
for 0(25%) of all signal events. Hence, we expect a fairly flat negative correlation between 
the two couplings. In the gluon-fusion production channel (and to some degree in the 
weak-boson-fusion channel) we can observe both decays to heavy gauge bosons, which 
means that the ratio of the two couplings protect the ratio of the two branching ratios. 
These two effects together with the general positive correlation described above again pull 
the gwwH — gzzH correlation into opposite directions and largely cancel. This is confirmed 
by Fig. P where the quality of the gbbH extraction is reduced and the positive correlation 
becomes dominant. 

In the third panel we show the correlation of the Yukawa couplings gbbH and grrH, J-e. 
the ratio of the two down-type Yukawas. In many models, like in type-II 2HDM this ratio 
will be invariant under leading modifications for example by an additional factor tan/3. 
For LHC measurements, these two Yukawas are not directly linked, if we neglect the two 
numerically marginal contributions to g-y-yH- However, we see the over-all scaling of all 
couplings also between these two. 

Last, but not least we show the correlation between the top Yukawa coupling and the 
on-shell top mass. The latter we can measure at the LHC to 0(1 GeV) and is already 



measured at the Tevatron to 1.2 GeV fT^]. We see that at this level a possible correlation 
of the top Yukawa coupling in the numerator of the effective gluon-Higgs coupling and the 
loop suppression by the top mass plays hardly any role. 

Based on Figs. |l] and ^ we can already conclude that for this kind of analysis Bayesian 
probabilities are less useful to understand the correlations around the best-fit parameter 
points and that alternative global maxima will likely not play a big role in our analysis. 
Therefore, for the rest of this paper we limit ourselves to (cooling) profile likelihoods. 
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Figure 3: Profile lilcelihoods; the three rows correspond to (1) only a ggH effective coupling, 
(2) only a "/jH effective coupling, and (3) both of them. The corresponding results without any 
effective coupling are shown in Fig. [l]. 



In Fig. ^ we study in detail the impact of the additional effective QggH and g-^yH ■ Both 
of them appear in the Standard Model, but they are completely described by the tree-level 
QttH and gwWH couplings, with minor contributions from the bottom quark loop. Physics 
beyond the Standard Model will always appear in these loops, in weakly interacting as well 
as in strongly interacting models. 

We start from the discussion of Fig. || which indicates that without effective couplings 
there is a clear correlation between guH and gwwH, which are the two loops contributing 
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to the one-loop H^^ coupling. The question is what happens if we allow for an additional 
effective coupling as shown in eq.( |3.2[ ). In the Standard Model, such small couplings exist, 
and in a general TeV-scale physics model they can become strong. In Fig. ^ we first add an 
additional gggH- This is at the expense of one of the measurements of the top Yukawa, since 
the dimension-five ggH operator in the gluon-fusion production cross section is the main 
source of the guH measurement. This is confirmed by their correlation or anti-correlation, 
depending on the sign of the top Yukawa. Both signs of guH are in principle allowed, 
but their degeneracy is broken by the ^^H coupling. Since the sign of the sub-leading 
bottom Yukawa couplings is poorly constrained this implies that little changes as long as 
the top Yukawa couplings is still determined by the g-yyH- The correlation between gggH 
and grrH we show as an example how these correlations are fed through the parameter 
space, combining the effects of loop contributions and production and decay rates. 

Adding only an effective g'yyH immediately allows for an equally likely negative top 
Yukawa coupling. Again, we see a weak correlation between g^yyH and the top Yukawa, but 
since the latter is well determined through gggH it is quantitatively less relevant. We also 
see a correlation between g^yH and grrH, with four equally likely branches representing 
the two possible signs of gun combined with the modulo-two degeneracy of A^^h = —2 or 
zero. 

Finally, we allow for both effective couplings. The correlation between both effective 
couplings now shows eight distinct maxima, with modulo-two steps in AggH and A^^h 
and a sign switch of the top Yukawa coupling. The left pair of solutions in the gggH~ 
g-yyE correlation has a (correct) positive top Yukawa coupling and a shifted AggH ~ —2. 
The central two pairs include the Standard Model solution with positive top Yukawa and 
^ggH ~ as well as a negative top Yukawa coupling combined with AggH ~ 0. In the right 
pair the negative top Yukawa coupling is compensated by AggH ~ 2. The same features we 
see in the one-dimensional likelihood for the effective Higgs-gluon coupling. Note that these 
solutions are only degenerate because the general effective couplings remove any sensitivity 
on the signs of the different Higgs couplings — confirming that the sign-flip symmetry is 
only broken by interferences in the loop couplings. 

4.2 Theory errors 

In general, to compare different hypotheses with data the central values of a parameter fit 
are useless without appropriate error bars. Therefore, also in this Higgs sector analysis we 
have to ensure that the error treatment is correct and that we understand the impact of 
different errors. 

While the proper treatment of statistical and (correlated) systematic experimental er- 
rors is well established and standard in all analyses, the treatment of theory errors deserves 
a closer look. In Fig. ^ we compare the impact of theory errors on the (hypothetical) most 
constraining scenario we can imagine: a parameter set without effective couplings and 
assuming a full 300 fb~^ of integrated luminosity. Obviously, we expect theory errors to 
have a sizable effect on the results. The second scenario for which we study the effect 
of theory errors includes both effective couplings and a lower luminosity of 30 fb~^. As 
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Figure 4: Profile likelihoods without (left) and including (right) theory errors. The top rows 
assume 300 fb~^ and no effective couplings, the bottom rows 30 fb~^ and both effective ggH and 
jjH couplings. 



discussed in Section 4.1, including the two additional effective couplings has a strong effect 
on correlations and on alternative minima. 

Including theory errors on the right-hand side we clearly see the flat region with max- 
imum log-likelihood around the correct central value spanning a few bins. This is a direct 
consequence of the Rfit scheme. For each measurement there is a flat region around the 
central value with the size of the corresponding theory error and likelihood one. This 
multi-dimensional box translates into a region in parameter space, i.e. we can slightly vary 
each parameter around its central value without any penalty in the likelihood. Outside this 
flat region the curves should not drop equally fast in both cases. This is indeed visible, in 
particular for the 30 fb^^ case where the fall-off is not as steep as for 300 fb^^. For example 
the tail for large grrH looks exactly the same. It is interesting to note that theory errors 
also increase the likelihood of the alternative solution Auh = — 2 for forbidden effective 
couplings. This alters the effective coupling of the Higgs to photons which within theory 
errors we can slightly shift without increasing x^- 

Theory errors can affect background predictions in addition to signal rates. However, 
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Figure 5: Profile lilcelihoods for tlie standard analysis (top row, copied from Fig. an increased 
ccH coupling with no additional contribution to the Higgs width (second row), an increased ccH 
coupling and a scaling factor for the total width (third row), and finally this scaling factor only 
(single bottom panel). All experimental and theory errors included for low- luminosity running. 



most of the LHC Higgs searches rely on control regions or side bins. The classical example 
is the two-sided side bins in the — > 77 search where the theory error does not have 
to account for higher-order corrections of the background rate. Therefore, we neglect the 
effect of background theory errors. The systematic errors on background rates are listed 
in Table |. 

4.3 Unobserved or invisible? 

There are two kinds of Higgs decays which we have not yet considered in this analysis. First, 
the Higgs boson can decay into invisible particles, like dark-matter agents, which would 
appear as missing transverse energy in the detector. As discussed in the introduction we 
skip this kind of analysis. 

What is more interesting is an actually unobservable Higgs decay, i.e. a Higgs decay 
which at the LHC we could not see. Looking at the particle content of the Standard 
Model, any decay into light leptons would be immediately visible. Shifted couplings to 
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third-generation fermions are already part of our parameter set, as is the couphng to 
gluons and the electroweak gauge bosons. One thing that could still happen are Higgs 
decays to light quarks, which would completely vanish in the QCD continuum. 

As an example, we consider an increase of the charm Yukawa coupling to 15.4 GeV, 
which corresponds to a branching ratio of roughly 80%. To be consistent, we also include 
this coupling in the rate for inclusive Higgs production. Identifying the light-jet coupling 
with the charm Yukawa is conservative in the sense that it minimizes the production rate 
for a given light-jet partial width because of the small charm parton densities. 

In the parameter extraction, the first effect of this additional coupling is the increase 
of the inclusive production rate. This effect is stronger than suggested by just comparing 
QggH and QccH, because the latter induces a tree- level decay. This rate enhancement can 
be counter balanced by reducing the branching ratios, i.e. by reducing all relevant Higgs 
couplings to final-state particles. This implies very small rates from weak boson fusion, 
driving them deep into the Poisson region for 30 fb^^. The effect on the parameter ex- 
traction we see in Fig. ^: in the second row the most likely points of all couplings are 
shifted by roughly the same amount. In addition, the errors are increased significantly, in 
particular for grrH- The separate peaks for QggH, corresponding to the different signs of 
the individual couplings, get smeared out into a large band of possible solutions. 

The consequence of such a result should be a more careful analysis, now including 
a free Higgs width. Such a shift shown in eq.( |3.4D can be implemented as a decay to 
an unobservable new-physics particle which does not appear inside the proton. Known 



examples for such effects are unobservable four-body Higgs decays |71|: The result we 
show in the third row of Fig. ^. We always find a best-fitting parameter point and the 
error bars are hugely inflated, in particular to large couplings. For the Standard Model 
couplings only the region around A = —1, where the coupling vanishes, is excluded. The 
smallest possible solution we obtain when the additional contribution to the width is zero. 
In that case we reproduce the slope from the no-additional-contribution scenario, modulo 
statistics effects. The region of large couplings is similar to the scaling of a free total Higgs 
width discussed in Section ^ Any increase can be compensated by a corresponding increase 
of the invisible decay mode. The predictivity of such an analysis is therefore rather limited. 

There is a second scenario we can test at this stage: if there appears to be a problem 
with the total width we can fix all Standard Model couplings to their respective nominal 
values and only allow for an additional experimentally unobservable decay mode. This 
leads to a very simple one-dimensional analysis, where only Ar is varied. Technically, we 
identify this increase with a large ccH coupling, but without including its contribution 
to the Higgs production rate. Under this assumption we can determine the size of the 
additional contribution to the width, as shown in the bottom line of Fig. ^. 

Finally, there is a third way to establish an upper bound on additional contributions 
to the Higgs width. Following the previous subsection we cannot simply add an additional 
contribution to the total width, because this would eliminate the possibility of determining 
the absolute scale of the couplings. On the other hand, following the argument in Section 
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Figure 6: Profile likelihoods for the standard parameter analysis (top row, copied from Fig. |^), 
and an additional contribution to the total width while at the same time keeping gwwH fixed 
to its Standard model value (bottom rows). All experimental and theory errors included for low- 
luminosity running. 



we can trade the assumption on the Higgs width for a limit in any other coupling, for 
example fixing gwwH to its Standard Model value. As we will see in Section this is 
the best-determined coupling of the parameter set. In models with more than one doublet 
this means that we either assume a rapid decoupling of the additional heavy states or close 
enough Higgs modes that we cannot resolve them in the transverse mass reconstruction of 
the WW final state (7|. 

In Fig. I we see that this indeed retains well-defined solutions. The errors on the 
different parameters stay at a similar level as the original analysis, with gi,hfj and grrH 
being measured mildly more precisely. Both of them are determined in channels with a 
production-side gwwH- The impact of the Z boson replacing the W in these channels is 
always sub-leading. The additional coupling ranges around zero with a standard deviation 
of roughly twice the Standard-Model width. This error corresponds to the error in the 
bottom Yukawa coupling, which is the main observable partial width and well determined. 



The results for these different scenarios give us some confidence in the robustness of 
our weak-scale analysis. While the LHC cannot make any direct statement about the total 
Higgs width (as discussed in Section ^), effects from unobserved decays or shortcomings 
of our effective Lagrangian hypothesis will surface as unexplained effects in the parameter 
extraction. These effects then trigger for example modifications of our ansatz and can be 
studied. 



4.4 Coupling ratios 

Given the sizable correlations in the Higgs-sector parameter space we could try to define 
better-suited parameters than individual Higgs couplings. In a way, this is in analogy to 
the extraction of super symmetric masses from cascade decays, where we can extract mass 
differences much more precisely than the actual masses. In the Higgs sector the hope is that 
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Figure 7: Profile likeliiioods for the absolute couplings (top row, same input as in Fig H), for the 
logarithm of the ratio of the gzzH, gtm, ^ttH (center row, left to right), gbbH, 9ggH and g-f^H 
(bottom row, left to right) couplings and the gwwH coupling. All couplings/ratios are normalized 
to their Standard Model value. We include all experimental and theory errors for a running of 
30 fb-i. 



ratios of couplings are better suited for an extraction from LHC measurements, for example 



coming from different decays but the same production process (or vice versa) 11]. In 
this situation some of the systematic and theory errors will cancel. The same trick allows us 
to extract flavor-physics limits from data in the presence of large uncertainties in low-energy 
QCD predictions or form factors. 

In general, we might expect for example coupling ratios of 1/2 and 2 to be symmetric 
and equally likely. However, naively showing the coupling ratio would tilt this distribu- 
tion, so in Fig. |^ we instead show the logarithms of the absolute values of coupling ratios, 
normalized to their respective Standard Model values; gwwH serves as the general nor- 
malization. For the couplings present in the Standard Model we see the expected peak 
around zero. The secondary solution of lS.ttH = —2 also gets mapped to zero by our for- 
mulas for the ratios. For the effective couplings gggu and g^p/H we see secondary peaks, 
originating from the alternative solutions in Fig. |3[ For the ggH coupling both the peaks 
at lS.ggH = and —2, where the latter corresponds to a top coupling with the correct sign 
and over-compensation by the effective coupling, appear at zero. The peak at ^ggH = 2, 
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compensating for a negative top coupling, appears just above one for the logarithmic ra- 
tio. In the 77i? case we see in principle the same effects. As the top-quark contribution 
is sub-leading we now have four solutions. They appear as three different peaks in the 
logarithmic ratio, because again A^^// = and —2 coincide at zero deviation. 

Both effective couplings also show a softer shoulder below the peak than above. The 
same effect we see for gzzH and guH, albeit less pronounced. This is due to the fact that 
reducing the couplings gives us a well-defined worst case of zero coupling, i.e. no signal 
events. In contrast, increasing a coupling we can generate arbitrarily large numbers of 
events with an arbitrarily small inverse likelihood. 

To answer the crucial question, we show the coupling constants itself in the top row 
of Fig. 1^ For gzzH we see hardly any difference between the coupling and the coupling 
ratio. The main measurement of this coupling arises from the gluon-fusion initial state 
where W bosons do not play any role. For quh the situation is similar; both, guH and 
9WWH are determined from a multitude of measurements so there is no particular benefit 
in looking at this ratio. In contrast, for grrH we can determine the ratio more precisely 
than the absolute coupling, because of the two final states in weak-boson-fusion production. 
The most striking improvement we find for gbbH^ where we exploit the correlation between 
the different occurrences of gwwH and the total width, dominated by bottom decays. 
The detailed discussion of errors on the ratios we defer to Section |5]^; it will give us an 
improvement on the extraction of typically roughly ~ 10%. 

5. Higgs couplings 

Before we can start determining the errors on the Higgs-sector parameters given a complete 
simulated data set at the LHC, we have one final step to implement. Until now, we have 
always tested our model on a set of perfect measurements, i.e. these measurements sit on 
their central value as predicted by the Standard Model without any smearing. As long as 
we do not extract error bars from such a data set this assumption is safe. However, for 
the final word on errors we need to move the signal and background rates away from the 
central value simulating an expected statistical fluctuation. 

5.1 Smearing the data 

One- and two-dimensional proflle likelihoods for the Higgs parameter extracted from a 
smeared set of measurements we show in Fig. ^. This smearing of the one pseudo- 
measurement includes all experimental but no theory errors. We see that the behavior 
of our results does not change significantly whether we use the true or smeared data set. 
The best-fitting points are slightly shifted away from the correct Standard-Model values. 
Moreover, the peak structures get broader, as the different measurements now attempt to 
pull the parameters into different directions. This indicates that any kind of error estimate 
from these likelihood distributions has to be taken with a grain of salt and we better rely 
on toy experiments Q to extract reliable error bars for the model parameters. 
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Figure 8: Profile likelihoods for unsmeared (top, taken from Fig. ||) and smeared (bottom) mea- 
surements assuming 30 fb~^. We include both experimental and theory errors in our analysis. 

In Table ^ we list the five best-fitting parameter points, obtained from an additional 
Minuit run starting from the best points of our (weighted and cooling) Markov chains. The 
first point is the closest to the Standard Model, where all couplings have the correct sign. 
All alternative solutions appear after flipping the sign of one or several of the couplings 
via A = —2. Only in the j'jH effective coupling, where the W and the top interfere 
destructively does the sign play a role. In the second column the sign of the two effective 
couplings is changed. The additional contributions to the effective couplings just over- 



best solutions 





0.12 


0.11 


0.11 


0.09 


0.09 


^ZZH 


0.12 


0.13 


0.11 


0.10 


0.09 


AttH 


-0.29 


-0.29 


-0.30 


-1.70 


-1.70 


^bbH 


0.17 


0.16 


-2.15 


-2.12 


-2.11 


^ttH 


0.02 


0.03 


0.03 


-0.10 


-0.10 


AjjH 


-0.01 


-2.43 


-2.42 


-0.38 


-2.76 




0.47 


-1.85 


0.34 


-0.48 


-0.47 



Table 4: Best solutions for one (experimentally) smeared data set with 30 fb ^ . The Higgs mass 
is 120 GeV. Both experimental and theory errors included. 
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compensate the Standard-Model parts. In the third column the bottom Yukawa coupling 
has a negative sign and this is reflected in a slight decrease in the effective gluon coupling. 
When we flip the signs from both the top and bottom quark in column four the ggH 
contribution is back to its original absolute value. In addition, we see the shift in the 
effective photon coupling, which has to compensate for the now constructive interference 
of the W and t loops. The same shift appears in the fifth column where the effective photon 
coupling over-compensates the Standard-Model contributions. 

Note that none of these solutions are 'true' alternative solutions, i.e. minima in 
not correlated to the best-fit values and with a distinctly different underlying set of model 
parameters. This by hindsight justifies our choice of profile likelihoods for this analysis, as 
compared to Bayesian probabilities. Practically it also means that for the remaining part 
of this paper we focus on the 'most physical' parameter point, i.e. we ignore secondary 
points with Ajjh ~ — 2. 

5.2 Best fit with error bars 

In Table ^ we show the errors on the extraction of Higgs coupling parameters. These 
errors we obtain from 10000 toy experiments, each smeared around the true data point 
including all experimental and theory errors. Besides calculating a root-mean-square error, 
we histogram the best fits for each parameter and extract (Jsymm using a Gaussian fit. As we 
do not expect the errors to be symmetric, we also fit a combination of two Gaussians with 
the same maximum and the same value at the maximum, but different widths. Depending 
on whether we are below or above the maximum we use the first (o"neg) or the second 
Gaussian (upos). For AzzH we see that as we approach —1, or vanishing coupling, the 
histogram goes to a constant value and is not well-fitted by a Gaussian. Therefore, for the 
lower branch of this coupling we fit only the central part within one standard deviation. 
Using a different measure for the shape of the error distribution we also show the standard 
RMS values for all couplings. They are systematically larger, owed to individual toy 
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Table 5: Errors on the measurements from 10000 toy experiments. We quote errors for Standard 
Model couplings only and including effective ggH and ^jH couplings using 30 fb~^ of integrated 
luminosity. The different a measures are defined in the text. 
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experiments far from the best-fit points. Only for symmetric Gaussian behavior we can 
expect these three error measures to coincide. 

First, we see that a correlation between the mass measurements and the couplings plays 
hardly any role. The errors on the Higgs mass are symmetric and the values correspond 
to our input data. 

A bit surprisingly, gwwH is basically unchanged whether we allow effective couplings 
or not. This coupling is the main contribution to the effective g-y-yH, so we would expect 
that allowing for additional effective coupling should remove the decay to photons from 
the gwwH determination. The results indicate that the accuracy from the remaining 
measurements is sufficient to determine gwwH- 

For the ttH coupling the situation is different. Besides its small tree-level impact it 
is the dominating contribution to gggH and a sub-leading contribution to g-y-yH- Hence, 
allowing effective couplings increases the errors. With and without effective couplings a 
Gaussian does not describe the errors accurately. In fact, including effective couplings we 
see a flat part around the Standard Model value, originating from our flat theory errors, 
and then a slightly steeper exponential fall-off. 

The ttH and bbH couplings are both strongly linked to gwwH, because of their 
respective production modes. The tt- associated production channel has too large errors to 
play a significant role. Consequently, we do not see any dependence on the existence of 
effective couplings. 

The ZZH coupling shows a particular effect when we include the effective couplings; 
its error is decreased. Such a behavior is counter-intuitive, because naively we would 
expect that such an additional coupling can always be considered a nuisance parameter 
which we need to remove. No matter if we integrate it out or project it away, the resulting 
error should form some kind of envelope. However, the situation changes once we include 
correlations. Let us consider two couplings gi and g2 (with a strong positive correlation) 
contributing to the same observable as g'f^l- With a combined shift the relative errors 5i 
and 82 on the couplings yield a total shift fi'i^lC'^i ~^ ^2) on the observable, limited by the 
measurement. For example, if due to the positive correlation both terms in the parenthesis 
are strictly positive, the upper limit on the observable now splits into two contributions. 
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Table 6: Errors on the ratio of couplings for a 120 GeV Higgs, corresponding to Table ||. The 
coupling ratio is normalized to the Standard Model value according to eq. 5.1. For the Standard 
Model dataset used in this study the central value of ^jjH/wwH is zero, therefore only the errors 
are quoted. 



-28- 



one for each parameter. In contrast, if we fix 6gi = 5g2 alone accounts for the complete 
deviation in the observable. 

The main source for determining gzzH is gluon-fusion production with subsequent 
Higgs decay into two Z. Without any effective couplings only the top Yukawa coupling 
determines the production rate, which is constrained by many other measurements. To 
account for the measured event rate we need to adjust gzzH accordingly. With additional 
contributions to gggH the top Yukawa coupling becomes largely irrelevant and the produc- 
tion and decay couplings of the gg ^ H ^ ZZ channel form exactly the combination 



discussed above. From Section LI we know that they are positively correlated, so indeed 
the error on /S.zzH can decrease. 

In Table ^ we show the errors on ^jjuiy/y^H, defined in the following equation as the 
deviation from 1 of the ratios of the coupling constants to the WWH coupling normalized 
to the Standard Model value: 

/ N SM 

— {^^) {l + ^nH/wwH) (5.1) 
gwwH \9wwhJ 

A non-zero central value for IS.jjjjiT^-^^ would arise in new physics scenarios. Thus 
A 

jjH/wwH provides information on a possible shift of the coupling ratios as well as the 
error on the coupling ratio. The central value of lS.jjH/wwH is zero for the Standard Model 
dataset used in this study, therefore we discuss only the errors on IS.jjHj\Yy/H. 

Table m confirms the qualitative results from the profile likelihoods in Section 44. The 



error on gzzH even increases once we form the ratio. The main determination mode is 
via the gluon-fusion initial state, which is independent of gwwH- An additional constraint 
enters via the subleading contribution to weak boson fusion, where gwwH and gzzH occur 
additively, so forming the ratio exacerbates deviations instead of decreasing them. For 
the top-quark coupling the effect of the effective photon coupling becomes clearly visible. 
Without this additional parameter guH and gwwH are linked and consequently the error 
on the ratio is significantly smaller than the error on the coupling itself. The additional 
g^^H breaks this correlation and we no longer gain anything by using the ratio. For all 
other couplings we observe more or less significantly smaller errors for ratios of couplings, 
because the production-side gwwH enters the determination of the decay-side couplings. 
Only the correlation between the well- measured gwwH and ghbH sticks out; it is due to the 
appearance of the total width in all rate predictions, which leads to strong correlations as 



discussed in Section 4.1. 



A note of caution should be added: the error on the error determination of the ratios 
of the coupling constants is non-zero. This is due to the fact that the coupling ratios have 
been determined from the coupling central values for a given toy experiment. Then for each 
toy experiment the ratio of the couplings calculated and the distribution of these ratios 
analyzed. The flat theory errors do not ensure a cancellation of correlated theory errors in 
these ratios. As the errors for the integrated luminosity of 30 fb~^ are dominated by the 
expected experimental error (78% to 94% of the total error as shown in Table [l|), we defer 
the development a more sophisticated treatment ensuring the full cancellation to a later 
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Figure 9: Profile likelihoods including the WH/ZH channel (upper row), WH/ZH with reduced 
sensitivity (center) and only ttH (lower) with a decay H —f bb. All experimental and theory errors 
included for low- luminosity running (true data points). 

date. We have checked however that by setting the theory errors to zero (implying a full 
cancellation of all errors) the error on the coupling ratios is improved by at most 10%. 

5.3 Bottom Yukawa coupling 

As mentioned in Section ^ the measurement of the bottom Yukawa coupling is crucial 
to the success of this parameter analysis; without it we would hugely underestimate the 
total width of the Higgs boson and scale down all couphngs correspondingly to obtain the 
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measured production and decay rates. For a long time the LHC experiments expected the 
ttH production process to be best suited for Higgs decaying to bottoms. Such a complex 
signature is more promising than simply looking for inclusive bb production or weak-boson- 
fusion production, because it includes more information to reduce backgrounds. On the 
other hand, several bottom jets lead to combinatorial errors, which turn out to largely kill 
this signature. Recently, there has been progress in understanding the structure of bottom 
jets coming from a Higgs resonance as compared to continuum QCD production. These 
two sources can be distinguished using information from the merging of subjets inside the 
jet algorithm used. Triggering requires for example additional leptons in the final state, 
which makes VH the preferred production channel. While these results have not yet been 
confirmed by detailed simulations, we nevertheless use them to illustrate the impact of the 
H ^ bb channels on our analysis. 

In Fig. P we first show results using the published WH, H ^ bb results at face value. 
This is the signature we base our entire analysis on. To be consistent with the other 
channels, we use the rates as in Ref. but included our own systematic and theory 
errors. 

In the second line we show results with a 50% reduced signal rate. As expected, all error 
bars are slightly increased, but the features of the distributions are virtually unchanged. 
Focusing on we see that even with the reduced sensitivity a precise measurement is 
possible. The slight degradation of this coupling (and with it of the total width) then 
propagates into all remaining couplings. 

Once we remove the WH/ZH sub jet channel from our analysis a determination of 
the bottom Yukawa coupling becomes impossible. Its error bar extends to a continuum 
0(150%) around the input value, and only a zero coupling is significantly less likely. The 
one measurement to determine g^bH is top-quark-associated Higgs production. However, 
for this channel a measurement of zero events only requires a fluctuation of 0.7 standard 
deviations. An additional constraint enters via the total width, but we already know that 
there is a scaling symmetry. And finally, for the bottom contributions to the effective 
QggH and Q'yyH to play any role we need an enhancement of the bottom-Higgs coupling far 
beyond those considered here. 

Aside from a broader peak structure of Abf,H we also observe considerably higher 
shoulders towards large values of the bottom Yukawa coupling. In this regime the branching 
ratio into bottom quarks gradually approaches one, so any increase in the coupling hardly 
shifts the number of H ^ bb events. In the center of the AhbH distribution the Ax^ 
value for vanishing coupling ranges around 1.5 in the absence of the sub jet analysis. The 
ttH channel cannot give such a large value, which means there is an additional indirect 
contribution from the total width. 

Precisely due to the increased uncertainty on the total width the otherwise well deter- 
mined gwwH suffers visibly. Not surprisingly, this increased error bar immediately propa- 
gates into QttH because of the correlated production process. In contrast, the top Yukawa 
coupling is virtually unchanged: first, the ttH production channel of bottom quarks con- 
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tributes nothing to the top Yukawa measurement, and secondly the bottom contribution 
to gggH is smah enough to still vanish in the error. The same is true for the effective QggH 
shown in Fig. ^. Both these couplings appear only on the production side of our measure- 
ments, so the bigger uncertainty on g^hH will already be absorbed by the respective decay 
channel. The situation is different for the effective Q'yyH, where the broader range of the 
^WWH requires a larger variation of A^^jj to fulfill the constraints in particular from the 
well-measured inclusive ^ 77 channel. 



As already alluded to in Section 4T, with most measurements significantly degraded 
we now observe the proper correlation between gwwH, gzzH and grrH from the weak- 
boson-fusion channels. Moreover, we see the positive correlation from the two observed 
decays in gluon fusion. 

In Table | we show errors on the parameters for the three different scenarios. Again 
we fit only the central peak around the Standard Model solution. These numbers confirm 
the findings from Fig. |^ A very prominent shift we see of course in ghbH^ as this is directly 
affected by changes in the subjet analysis. There we also see a large fraction of solutions 
where the coupling vanishes. Other couplings are extracted with a reduced precision of 
50 — 100%. Only the error on gun is largely unchanged with remaining differences in part 
due to statistical fluctuations. 

5.4 Minijet veto 

One of the crucial ingredients to the weak-boson fusion production analyses is the central 



mini-jet veto |29| . The two tagging jets are emitted into the forward regions of the detector, 



while the central region stays free of jets. This behavior is reflected in the small next-to- 



leading order QCD corrections to this process |^, 77 1. Typical backgrounds as well as 
gluon-fusion production involving two jets do not share this feature. How well Monte- 
Carlo simulations describe the mini-jet veto probability is not yet established. Without 
any claim of a useful range we increase the uncertainty on the mini-jet veto probability 
from 5% to clearly unrealistic 100%. 

In Table ^ we see that even with this huge theory error the numbers hardly change 
at all. The fit of a symmetric Gaussian to the toy experiments consistently enlarges the 
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Table 7: Errors on the measurements including the WH/ZH channel (left, from Table||), WH/ZH 
with reduced sensitivity (center) and only ttH,H — > bb (right). We assume 30 fb"-'^ of integrated 
luminosity. 
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Figure 10: Profile likelihood for /^.ttH including the standard error of 5% (left) and an increased 
error of 100% (right) on the central mini-jet veto. We show results for true data points and 30 fb~^. 

errors, although the effect is mostly minimal. The varying behavior for the asymmetric 
errors can be traced back to statistical fluctuations. Only for grrH we see a significant 
change compared to our normal analysis. This coupling is determined exclusively via 
weak-boson-fusion channels, so an increase in the error directly affects this coupling. In 
Fig. |l^ we see a clearly finite likelihood for vanishing coupling. In the complete scan we 
even see a slight peak in this place, because due to the large errors it is now possible that 
the rates for the weak-boson-fusion channels fluctuate below the background. To illustrate 
a non-negligible effect of the increased error: we see that while the zero-signal solution is 
clearly excluded in the standard case, we now cannot even reach a 95% confidence level 
assuming an integrated luminosity 30 fb~^. 

While this discussion seems to indicate that the exact knowledge of the minijet veto 
probabilities is not crucial to our analysis, there are two caveats: first of all, we only 
study the effect of the minijet veto probability for the signal, while we assume that it 
can be measured precisely for the backgrounds. This assumption might or might not be 
valid once we face real data. In particular in tough analyses like invisible Higgs searches 
(which we do not use in this analysis) this assumption needs to be reconsidered. Secondly, 
we limit ourselves to low-luminosity running, which hurts the statistical pull of the weak- 
boson-fusion channels. The moment we consider for example 100 fb~^ of data this situation 
will likely change. 
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Table 8: Errors assuming the 5% standard error on the mini-jet veto (left, from Table ||), and 
assuming an increased error of 100% (right), using 30 fb^^ of integrated luminosity. 
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6. Beyond the Standard Model 



The main aim of the kind of Higgs sector analysis presented in this paper is to probe the 
nature of the Higgs sector in case that we only see a light scalar Higgs boson at the LHC. 
In this section we study the application of the Higgs parameter analysis and the ansatz 
described in Section ^ on a supersymmetric Higgs sector. We restrict the analysis to the 
standard without additional observables |65]. 

This ultraviolet completion of the Standard Model has the advantage or disadvantage 
that it contains a the decoupling limit in which the LHC might quite realistically be left 
with a single light Higgs boson with very similar properties to its Standard-Model counter 
part. Such an outcome is predicted for fairly small values of tan/3 < 10, where the Yukawa 
couplings of the heavy Higgs states are not sufficiently enhanced to warrant a direct obser- 
vation. In this regime the only hope to see a heavy Higgs at the LHC might be an on-shell 



decay H ^ hh with a subsequent decay for example to 6677 [37]. On the other hand, for 
small tan (3 the decoupling limit is approached far less rapidly, so that a careful analysis of 
all couplings of an observed light Higgs state might reveal significant deviations from their 
Standard-Model values. 

In this section we study two different scenarios: first, we choose a fairly low mass for 
the heavy CP-odd Higgs boson niA, so that we are not too far in the decoupling regime. 
In a second part we consider a gluophobic Standard-Model like Higgs, i.e. the top loop in 
the effective ggH coupling is largely canceled by a corresponding stop loop |7^. Just like 
in Section ^ we use true data points to disentangle the effects of the models from smearing. 
In addition, we use the same effective theory as described in Section ^ without including 
new Higgs search channels in the decays of new states and without including loop effects 
from particles seen elsewhere in Atlas or CMS. 

6.1 Supersymmetric Higgs 

As a first step of this application we need to replace all Standard Model measurements 



with the respective rates for a supersymmetric parameter point based on SPSla |75]. After 



evolving the original SPSla parameters from the high scale to the low scale using SuS- 



pect [76|, we change tan/3 to 7, At to —1100 GeV and ttia to 151 GeV. For this parameter 
choice gwwH and gzzH are significantly smaller than in the Standard Model. It also pre- 
dicts a light CP-even Higgs mass of 120 GeV. The true values for all Higgs couplings we 
give in the first line of Table ^. Comparing these deviations with the Standard Model error 
bands in Table we can guess that only (75(,H and grrH will be useful individually, provided 





^WWH 


^ZZH 




AbbH 


^ttH 




^ggH 


rriH 


true 


-0.13 


-0.13 


-0.19 


3.27 


3.29 


0.19 


-0.28 


120.0 


errors 


±0.45 


±0.61 


±0.63 


±2.34 


±3.35 


±0.99 


±1.12 


±0.29 




-0.43 


-0.99 


-0.60 


-3.68 


-3.23 


-0.70 


-0.69 


-0.29 




+0.48 


+0.52 


+0.65 


+ 1.52 


+3.58 


+ 1.30 


+1.46 


+0.30 



Table 9: Couplings for the SPSla-inspired scenario. We give the true values of our input and the 
error bars assuming 10000 toy experiments and 30 fb~^ of data. 
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Figure 11: Profile likelihoods for the parameter point SPSla with light itla and modified At and 
tan /3 for 30 fb~^. We include both experimental and theory errors in our analysis. 



we can extract them in this scenario. While the couplings to gauge bosons are different in 
our supersymmetric scenario, their shifts with respect to the Standard Model values are 
all within their respective errors. 

For the fairly small value of tan/3 the heavy Higgs bosons H, A and will not be 
observed at the LHC, and certainly not using 30 fb~^ of data. The question is: can we 
distinguish such a scenario from a Standard-Model Higgs sector just looking for a light 
Higgs boson [73|; 



? 



Strictly speaking, we should account for the extended particle content of supersym- 
metry with an increased theory error on the perturbative predictions of the production 
and decay rates. However, in most production channels the quantum effects of new states 
are mass suppressed and hence within the errors quoted in Table. ^ With real 

LHC data this procedure would be an iterative process: first we would try to establish a 
deviation assuming the Standard Model particle content as the correct hypothesis. The 
general new-physics mass scale accounting for deviations from the Standard Model cou- 
plings will then feed back into our new physics analysis. In turn, taking information from 
other sectors into account we can refine the Higgs sector analysis. 

In the first line of Table ^ we can guess that essentially the four measurements with 
T leptons or bottom quarks in the final state will contribute to our likelihood. All other 
measurements have too low production rates, even consistent with a null result with 30 fb~^. 
We nevertheless include all measurements in our analysis for consistency. For the dominant 
measurements guH 9WWH and gzzH appear on the production side. The latter is a sub- 
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leading contribution to both weak boson fusion and the subjet analysis. For these four 
parameters we have four high-quality measurements, so we can hope to obtain a unique 
solution in this subspace when considering the complete set of measurements. 

In the one-dimensional distributions in Fig. |ll| we immediately see that the Standard 
Model hypothesis is not correct. For g^hn we do not see an upper limit on the coupling 
strength. This is a consequence of the branching ratio into bottom quarks which event at 
unity would be compatible with the data. This limit removes any kind of rate measurement 
involving H ^ bb from our measurements. The question remains if this leads to the free 
scaling of the total width with all couplings, as we have seen before. 

In the subjet analysis channel we cannot make the branching ratio into bottom quarks 
larger than one, which induces a lower limit on gwwH- On the other hand, we cannot 
decrease the branching ratio significantly, because this would require a large value oi gwwH- 
Weak-boson-fusion production with subsequent decays into a pair of W bosons prohibits 
this solution, giving us an upper limit on gwWH and inducing the lower limit on gbbH in 
Fig. 0. 

The same argument we can now apply to guH and tt-associated Higgs production. 
The lower limit on g^bH induces an upper limit on guH, which we also see in our results. 
Additional information we obtain from the non-observation of other channels with top- 
quark initial states. As this measurement is much less precise than the subjet analysis, the 
allowed values of the top-quark coupling extend further to positive values. The lower limit 
again reflects the fact that the bottom branching ratio cannot be larger than one. 

The fourth parameter A-rrH is determined via weak boson fusion. The range in gwwH 
implies an allowed range of branching ratios into a r pair, and therefore of the ratio 
grrH/gbbH via the total width. Accordingly we see a strong positive correlation between 
these two couplings. This way, grrH can extend to arbitrarily large values, while we have 
a definite lower limit and the Standard-Model value is clearly excluded. 

For the remaining couplings we can extract limits from a non-observation in Higgs 
search channels. In g-y-yH we again see a scaling with the total width and therefore a strong 
correlation with gbbH- The boundaries of A^^h are simply given by our (arbitrary) limits 
on gbbH in this scan, as we see in the correlation plot. The excluded region around —1 is a 
normalization effect: the WWH and ttH couplings contributing to the effective coupling 
are smaller than in the Standard Model, while our additional contribution is normalized to 
the Standard Model value. Setting A^^h = — 1 would not give a vanishing total coupling 
but something too large. The bands we see are the regions where the additional contribution 
cancels the loop-induced coupling. 

For the gggH coupling we do not see a scaling with gbbH, because any effect of varying 
this coupling will already be compensated by the associated decay coupling. Including 
the tan/3 enhancement the bottom loop can now become large enough to give a relevant 
contribution to the effective coupling, so we then need to dial the additional parameter 
accordingly to cancel the combinations from both heavy quarks. Therefore, we do not see 
a sharp peak in AggH but a wide band of best-fitting points. 
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Given such a scenario where we see clear deviations from the Standard Model, the 
question arises at which level can we exclude it. Obviously, in the presence of correlations 
and non-Gaussian errors we cannot simply add the individual log-likelihoods. Instead, we 
determine the exclusion limits using the log-likelihood q as an estimator. 

The first question is how well we can rule out the Standard Model prediction given 
a data set dsusY consistent with the superysmmetric Higgs model predictions msusv- 
The purely Standard Model log-likelihood distribution q{dsM\'m'SM) usually gives us the 
desired confidence level as the integral over q from minus infinity to the central value of 
the msusY- However, we also have to take into account that the new-physics data can 
fluctuate within its experimental and theoretical errors. Basing the statement on one 
single supcrsymmctric toy experiment would be incomplete. Therefore, we first compute 
a 90% confidence level for the Standard Model hypothesis given Standard Model data q. 
In a second step we evaluate the log-likelihood distribution for data consistent with the 
super symmetric scenario, but still based on the Standard Model couplings g'(dsusYl'^SM)- 
The percentage of toy experiments giving g((isuSY|"T^SM) < q we find to be 77%. This 
means that of all our toy experiments assuming the supersymmetric Higgs sector described 
in Table IX, 77% are not described by the Standard Model within the 90% confidence level 
and can be ruled out accordingly. 

The second question is whether the new physics model is a better description of our 
new physics data (dsusY) than the Standard Model. The relevant information is the 
distribution of the difference of the two log-likelihood values for each toy experiment: Aq = 
^(c^SUSyI'tisusy) — O'C^^SUSyI'tism)- Note that for most cases we expect this difference to be 
positive. Because of the non-standard distribution of our log-likelihood the usual definition 
corresponding to Ax^ > (with s a number of standard deviations) does not hold. A 
detailed discussion on testing different hypotheses and confidence levels we postpone to a 
later paper, while at this stage we simply request Aq > qQ where qo is chosen of the order 
of the quality of the consistently supersymmetric fit. Of all our toy experiments we find 
that only 4% are better described by the supersymmetric hypothesis than by the Standard 
Model given this naive choice. 

6.2 Gluophobic Higgs 

The determination of the Higgs mass plays an important role in the 120 — 140 GeV mass 
range because the different branching ratios are steep functions of the mass. For a low- 
mass Higgs the best mass measurement comes from the decay into two photons, which 
is usually extracted from overwhelming backgrounds using two side bands. The excellent 
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^ggH 
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true 


-0.01 


0.00 


0.00 


0.14 


0.13 


0.28 


-0.76 


112.36 


errors 


±0.37 


±0.53 


±0.50 


±0.67 


±0.62 


±0.68 


±0.69 


±6.83 




-0.41 


-1.21 


-0.59 


-0.39 


-0.24 


-0.48 


-0.59 


-1.84 




+0.32 


+0.29 


+0.42 


+1.10 


+0.99 


+0.86 


+0.79 


+ 10.21 



Table 10: Couplings for the gluophobic Higgs scenario. We give the true values of our input and 
the error bars assuming 10000 toy experiments and 30 fb~^ of data. 
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invariant mass resolution of the detectors then allows for a Higgs mass determination with 
a dominant experimental error around 0(100 MeV). However, in models where the loop- 
induced Higgs couplings to gluons is strongly reduced this mass measurement gets degraded 
by lower statistics. A super symmetric Higgs can be a good example for such a possible 
challenge. 

Naively, we would first expect problems with the effective coupling to photons including 
superpartners. Aside from the usual W and top loop, charged Higgs bosons, sfermions and 
charginos contribute to the effective coupling of the photons to the Higgs boson. However, 
the effective g-^yH is unlikely to decrease significantly after taking into account current mass 
limits. 

This is different for the inclusive Higgs production rate from gluon fusion. Here, the 
top loop naturally cancels with the stop loop due to the opposite spin of the two states. 
Moreover, the current experimental bounds on the stop mass is weak enough to allow for 
such a scenario, commonly called a gluophobic Higgs. The effect of this supersymmetric 
scenario on the Higgs parameters we show in the first line of Table. Similarly to the 
supersymmetric model shown in Table ^ individual deviations from the Standard Model 
couplings will not be conclusive, and the sizable ^ggH will turn out to be dangerous, 
because it reduces the number of Higgs bosons produced altogether and with that the 
precision on all model parameters. 



The reduced rate for the inclusive i7 ^ 77 signal means that this channel is no longer 
suited for a precise determination of the Higgs mass. An alternative is a Higgs decay to 
ZZ, but for low Higgs masses and strongly reduced inclusive production rates it is even 
less promising. For an integrated luminosity of 30 fb~^ the weak-boson-fusion production 
process with a subsequent decay to photons predicts 13.5 signal events with a background 
uncertainty of close to 20 events, which means it will only contribute to a Higgs mass 
measurement at higher luminosities. The last method for determining the Higgs mass is 
weak-boson-fusion production with a decay into a r pair. The rr invariant mass can be 
reconstructed in the collinear limit, with a mass resolution around 8 GeV |47, 46], which 
we assume as the error on the Higgs mass in this scenario. 



In Fig. 12 we show a set of profile likelihoods. The parameter determination suffers 
from the effective absence of all gluon-fusion production channels and the larger uncer- 
tainty in the Higgs boson mass. On the other hand, in contrast to the supersymmetric 
scenario discussed in Section |6.l| we do not observe any severe holes in the analysis. All 
parameters which do not depend, directly or indirectly, on gggH are determined well, with 
an increased error in part from the Higgs mass measurement. This includes gwwH, gzzH, 
gttH, and grrH- Even for the AggH itself the Standard Model value of zero is excluded. 
The four solutions for correct and flipped sign of the top-quark coupling and zero or dou- 
ble contribution from the additional effective coupling are now grouped in two pairs, each 
forming a broad peak. This we also observe in the two-dimensional likelihood. 

The bbH coupling depends on gggH via the total width. In the Standard Model loop- 
induced decays of the Higgs into gluons have a sizable branching ratio. A wide variation in 
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Figure 12: Profile likelihoods for a gluophobic Higgs scenario at 30 fb ^. We include both exper- 
imental and theory errors in our analysis. 



^ggH significantly shifts the total width and needs to be compensated by adjusting QbhH-, 
so the decays into bottoms account for the correct rates. The production-side QggH directly 
impacts the measurement of Q'yyH- Both, the large range for the gluon coupling and the 
increased statistical error on the inclusive photon channel lead to a poorer determination 
of this coupling. 

Again, we compute the significance with which we can exclude the Standard Model 
hypothesis. Follow the prescription outlined in the previous section we find that for 46% 
of our toy experiments we can rule out the Standard Model at a 90% confidence level. 
For data consistent with a gluophobic Higgs, the new-physics hypothesis gives a better 
description than the Standard Model in only 2% of the cases. 



7. Outlook 

In this paper we have studied quantitatively how well LHC measurements can determine 



the parameters in the effective weak-scale Higgs sector [ 1 1 1 . SFitter ^ uses Markov chains 
to map the multi-dimensional Higgs parameter space onto a set of measurements which 
we expect for a 120 GeV Higgs boson at the LHC. The resulting exclusive likelihood 
map allows us to show correlations and to determine the error bars for individual Higgs 
couplings. For an appropriate description of these errors we allow for general correlated 
statistical, systematic and theory errors. Since we do not find distinct alternative maxima 
in the exclusive likelihood map, we rely on profile likelihoods (as compared to Bayesian 
probabilities) to describe likelihood distributions in one or two dimensions. 
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We find that Higgs couplings can be extracted with typical errors around 20 — 40% 
using an integrated luminosity of 30 fb~^ and properly simulating all errors involved. The 
different parameters are strongly correlated, and sign ambiguities in particular in the pres- 
ence of dimension-five operators occur. One main correlation between all channels is due 
to the total width, which also means that a reliable measurement of the bottom Yukawa 
coupling for example using subjet analyses [41| is vital for our analysis. Coupling ratios 
instead of individual couplings can have somewhat reduced errors, also depending on the 
treatment of the total Higgs width. 



Unobservable (not 'invisible') Higgs decays at the LHC can include a coupling to light 
quarks which also affects the inclusive production cross section. We discuss strategies of 
dealing with different scenarios of this kind. While the LHC is clearly not going to measure 
the total Higgs width directly, unexpected effects not contributing to the observed decay 
channels will be visible in our fit. 

Given our results for the error on the couplings it is unlikely that we will be able to use a 
general Higgs-sector analysis to distinguish between different decoupling model hypotheses 
(like the MSSM vs Standard Model), but drastic modifications like low-m^ supersymmetry 
or a gluophobic Higgs boson will be clearly visible. We quantify the confidence levels of 
distinguishing the Standard Model hypothesis from the respective new physics hypotheses 
in a two-dimensional plane, including a possible fluctation of the data as well as a variation 
of the model predictions, as predicted by the complete error structure. 
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A. Combining log-likelihoods 

For each experimental channel in our analysis we need to combine different sources of 
errors. The number of events is given by a Poisson distribution which for non-integer 
event numbers we compute as P(m, d) = e~"^m'^ /T{d + 1), with m the predicted number 
of events and d the number of measured events. Systematic errors follow a Gaussian 
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Figure 13: Comparison of our approximate (dashed green) with the exactly convoluted (solid red) 
log-likelihood. We show the absolute and the relative difference in for 5 (left) and 50 (right) 
events with a Gaussian error of 10%. 

distribution, since they are measured using large background samples. We assume Higgs 
signals and background measurements of the systematics to be independent. Therefore, 
their combined probability is the convolution of the individual parts. For two Gaussians this 
can be done analytically, leading to another Gaussian with the errors added in quadrature. 
The convolution of a Gaussian and a Poisson distribution, on the other hand, can only be 
evaluated numerically. This is practically not feasible for each step in our Markov chain, 
so instead we use an approximate analytic formula. 

The basic quantity for SFitter is a generalized x^j related to the log-likelihood by = 
—2 log L. First, we transform the Poisson probability into a likelihood Lp. In analogy to the 
Gaussian case we normalize it to L{d = m) = 1. Since the peak of the Poisson distribution 
is close to d = m — 0.5 our likelihood would exceed unity for de [m — 1, m]. Therefore we 
fix it to unity in this interval. This way the maximum of a combined Gaussian and Poisson 
likelihood is well defined and unique at d = m. Compared to an exact calculation this 
method slightly overestimates the likelihood over the whole range, thereby erring on the 
side of caution. This leads us to 

Pirn d) 

-21ogLp = -2 log ' \ = -2[((i-m)logm + logr(m + 1) - logr(d+ 1)1 (A.l) 
F[m, m) 

which in the Gaussian limit becomes — 21ogLG' = {d~ riif' ja'^ with the combined error a. 

Next, we combine the two log-likelihoods into a single form, guided by two limiting 
cases: the combination of two Gaussians should return the exact formula for the combined 
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Gaussian. Secondly, if one of the log-likelihoods becomes very large, the result should 
approach the other log-likelihood. The form 



logL log Lp logLc 2 log ^("^'^) (d-m)2 

P[m, m) 

fulfills both conditions. To test this setup we can compare our result with a numerical 
convolution for typical cases as they appear in our data. In Fig. |l^ we show results for 5 
and 50 events. The Gaussian error is 10%, i.e. 0.5 and 5 respectively. In both cases the 
difference of our approximation to the exact result is small. Furthermore, we very slightly 
underestimate x^, and consequently slightly overestimate the errors, erring on the safe side. 

Note that our results can differ from the naive approach of just multiplying the indi- 
vidual likelihoods arising from different errors for the same channel. This approximation 
is only valid for different channels in the limit of vanishing correlations between these. 
We illustrate this for two errors in one channel: let us look at a measurement with two 
Gaussian distributions with a difference between predicted and measured values x and an 
error of 1 each. Adding the log-likelihoods gives us = 2x^/1^ = 2x^. To compute the 
standard deviation we set = \ and solve for x = l/\/2. In contrast, simply adding both 
errors in quadrature yields x = ^/2, which means adding the log- likelihoods underestimates 
the error. This tendency becomes more pronounced when the two errors are not equal. 



B. Cooling Markov chains 

The Markov Chain algorithm is optimized for creating Bayesian marginalized log-likelihood 
distributions. This is analogous to an integration, i.e. not only the height of a peak is 
relevant for its contribution to the final result, but also its volume in parameter space. 
Correspondingly, the Metropolis-Hastings algorithm involves two steps: first, it suggests 
the next point with a probability proportional to the parameter-space volume of a peak 
structure. Secondly, it accepts or rejects it based on the height of the peak. The total 
probability is a combination of the two. 

In a (frequentist) likelihood approach we are only interested in the height of the peak. 
Only the value at the tip determines the profile-likelihood value. Following the above 
argument sharp and high peaks are not very efficiently probed by Markov-Chain algorithms. 
Consequently, we amend the algorithm so that it better scans the vicinity of points with 
large function values. A well known approach which does this is simulated annealing |79|. 



Combining this with the Metropolis-Hastings algorithm jS^l we modify the accept or reject 
condition to 

/ suggested ^ ^100/(j-c) ^ ^-gj^-j 
/previous 

where /suggested and /previous are the function values of the newly suggested and the last 
point in the Markov chain, respectively. The complete Markov chain we divide into 100 
equally long partitions numbered by j. r is a random number between and 1, and c is 
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Figure 14: Comparison of profile likelihoods obtained from standard (left) and cooling (right) 
Markov chains. 



the cooling factor. We find that c ~ 10 gives well-converging results for this analysis. As 
in the original algorithm a better point is unconditionally accepted. 

In the beginning r is almost entirely mapped to values close to r^^^^^ ~ 0, i.e. we 
accept many points and quickly scan the parameter space. As j increases we turn via the 
standard algorithm into a region where we accept almost only better points. This way 
we carefully scan the vicinity of the peaks and we obtain reliable values for the profile- 
likelihood distribution. If there are several peaks of comparable quality, different Markov 
chains will focus on different peaks. We are guaranteed to see these alternative maxima 
because we combine 30 individual Markov chains for the final result. 

To show the power of this modification we show the one-dimensional profile likelihood 



for the WWH coupling in Fig. |14[ It includes 30 Markov chains with 200000 points 
each. We see how the coarse standard algorithm places the peaks at a completely different 
location while the cooling Markov chain gets both the position and the height of the 
maximum right (compared to the input). 
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